Triangular anisotropies in Driven Diffusive Systems: reconciliation of Up and Down 
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Deterministic coarse-grained descriptions of driven diffusive systems (DDS) have been hampered 
by apparent inconsistencies with kinetic Ising models of DDS. In the evolution towards the driven 
steady-state, "triangular" anisotropies in the two systems point in opposite directions with respect 
to the drive field. We show that this is non-universal behavior in the sense that the triangular 
anisotropy "flips" with local modifications of the Ising interactions. The sign and magnitude of the 
triangular anisotropy also vary with temperature. We have also flipped the anisotropy of coarse- 
grained models, though not yet at the latest stages of evolution. Our results illustrate the comparison 
of deterministic coarse-grained and stochastic Ising DDS studies to identify universal phenomena 
in driven systems. Coarse-grained systems are particularly attractive in terms of analysis and 
computational efficiency. 



Driven steady-state systems are common in many fields 
of physics including device physics and materials process- 
ing, and are also found in many biological processes. The 
simplest example of a driven system is one in which a uni- 
form external drive is applied. For closed boundary con- 
ditions, a final equilibrium state is reached in which the 
external drive is balanced by other forces — as in a closed 
system in a gravitational field. For open boundaries, 
equilibrium will never be reached and a non-equilibrium 
steady state will continue "indefinitely" . Electrical cir- 
cuits provide prosaic examples of this. In both open and 
closed systems, the introduction of local interactions pro- 
vide a rich phenomenology that is only slowly being ex- 
plored. 

To model a driven system we must describe both the 
energetics and the dynamics. One of the simplest such 
model is a lattice driven diffusive system (DDS). This is 
simply an interacting lattice gas with biased motion in 
the direction of the uniform external field. Early work 
on this model concentrated on steady-state properties 
near the critical point, which survives from the under- 
lying zero-field Ising model. More recently there has 
been increasing emphasis on the approach to the steady 
state. This dynamical regime is independent of the open 
or closed boundary conditions, and hence is common to 
both. See the book by Schmittmann and Zia Q for an 
introduction to the literature. 

One of the most important questions we can ask about 
any model is whether the behavior that it displays is uni- 
versal. To whit: is the observed behavior seen in a broad 
class of systems, or is it specific to that precise model? 
For example, we can ask whether a lattice DDS and re- 
lated coarse-grained models display the same universal 
behavior. This question has proven remarkably contro- 
versial in DDS, and this paper aims towards a reconcili- 
ation between lattice and coarse-grained approaches. 

Near the critical point separating the low-temperature 
ordered and high-temperature disordered DDS phases, 
coarse-grained descriptions lead to analytic solutions of 



the steady-state structure j2j. Ising DDS simulations 
have been consistent with these results jj| on balance, 
though they have been performed mainly at infinite fields 
Q. [There have been claims that different approaches, 
though still coarse-grained, should apply in that infinite- 
field limit [[3|(|.] Away from the critical point, the focus 
has been on the approach towards the final steady state 
and the question has been whether any coarse-grained 
model recovers the phenomenology of a stochastic Ising 
DDS. In the ordered steady-state of both the lattice gas 
and related continuum descriptions, the domain walls 
align with the field direction and fluctuations are sup- 
pressed to a remarkable degree E[. The main difference 
occurs in the approach to the steady-state, and is easi- 
est to visualize in a system with a marked minority of 
one of the ordered phases (so-called "off-critical" sys- 
tems) . There, separated droplets nucleate and diffusively 
coarsens [for the zero-field limit see 0]. As the drops 
grow, they elongate in the field direction — ultimately 
forming the stripes of the steady-state. This phenom- 
ena is seen both in the lattice DDS p| and in continuum 
models ||. However, the drop shapes are triangular and 
point in opposite directions with respect to the field in 
the Ising and coarse-grained models (see, however, JTo[). 
As seen in Fig. [j], the tips of the triangles point against 
the field direction (up) in the lattice model, while they 
point with the field (down) in the continuum simulations. 
The same triangular anisotropy, and discrepancy, is seen 
at other volume fractions as well. 

Does this indicate that either the Ising DDS mod- 
els, or the coarse-grained models (or both!) are non- 
universal? Either would be less than ideal, since stochas- 
tic (Ising) models are needed for temperature-dependent 
studies near the critical point, and deterministic coarse- 
grained models are both analytically tractable and more 
computationally efficient at lower temperatures. An at- 
tractive resolution would be that there are regimes of pa- 
rameter space in both stochastic Ising and deterministic 
coarse-grained models which exhibit either sign of trian- 



1 



TDGL 




- -• >.V 



50000 MCS 



t=1000 



I 



FIG. 1. Evolution of domains for an off-critical quench 
in a lattice gas with nearest neighbor hops (left) and for a 
TDGL coarse-grained model with isotropic mobility (right). 
The field E biases the motion of the particles (dark) down- 
wards. For the lattice gas the tips of the triangular domains 
point "up", opposite the field, for the coarse-grained model 
the domains point "down". 



gular anisotropy, down or up with respect to the field, 
with the previously studied models being in the two dif- 
ferent regimes. We take that as a working assumption, 
and try to flip the anisotropy by varying model param- 
eters in the Ising and coarse-grained approaches. We 
report partial success, with asymptotic reversal in the 
Ising DDS and reversal up to intermediate times in the 
coarse-grained DDS, and we have by no means exhausted 
the phase-space of model parameters. 

Even if we achieve full success, fundamental ques- 
tions remain. Is the triangular anisotropy simply a non- 
universal amplitude, or does it lead to different phe- 
nomenology? How do we understand the temperature 
and field dependence of the anisotropy, and how does 
the observed anisotropy in the dynamic correlations re- 
flect anisotropies due to the surface tension, particle mo- 
bility, and the external field? Reconciling the Ising and 
coarse-grained approaches is simply the start of the story. 



Lattice Gas Dynamics 

The basic discrete DDS model |],||§ is an extension 
of an Ising model with conserved Metropolis dynamics: 
particles hop with probability 



where the energy difference AH — A_ffi S i ng + Ai/pioid 
includes an applied field. Restricting ourselves to a two- 
dimensional square lattice, AHi s - lng is the standard near- 
est neighbor Ising Hamiltonian in which a particle inter- 
acts with its four nearest sites. [We absorb J/ks into the 
temperature, which we then measure with respect to the 
zero-field Ising T c ] AH-pi C \d is +E if the particle moves 
one lattice unit opposite the field direction and — E if 
the particle moves in the field direction, where E is the 
field strength. The dynamics are conserved, so particles 
hop rather than being created nor destroyed. Restrict- 
ing the hops to nearest neighboring sites, Alexander et 
al. H found that the domains formed upward pointing 
triangles, as shown in Fig. |[ 

Natural generalizations of this well studied DDS model 
include looking at different lattices [triangular, hexago- 
nal, and Kagome in 2d], rotating the field away from a 
lattice direction |TT| , allowing hops and/or interactions 
with further neighbor particles, and allowing anisotropic 
hop rates and interactions. Universal results on the basic 
model should be robust to these sorts of microscopic dif- 
ferences, as differences will certainly be entailed by exper- 
imental realizations. [These changes will affect variously 
the anisotropic particle mobility and interfacial surface 
tension of any coarse-grained representation.] In this pa- 
per we allow next-nearest-neighbor hops, where we treat 
all eight immediate neighbors with equal weight. We la- 
bel these "nnn" dynamics, in contrast to "nn" dynamics 
where hops are restricted to the four nearest-neighbors. 



Coarse-grained Dynamics 

The simplest coarse-grained dynamics is the time- 
dependent Ginzburg-Landau (TDGL) model with a field. 
The free energy is given as 



F[4>] = J 



dr 



f{4>(r)) + \W4>\ 2 + Ezdp 



(2) 



where 4>(r, t) is the order parameter and the field E points 
down toward lower z. Within a uniform phase, we use 
the following Flory-Huggins type free-energy density: 

}{4>) = (1 + cj>) ln(l + cf) + (1 - 4>) ln(l - 4>) ~ \<f- 

The system will phase separate for a > 2 with coexistence 
values depending on a. [For a k 2 this recovers a more 
familiar (j) 4 free energy] This choice of f((f>) forces \<f>\ < 
1, which simplifies the treatment of the particle mobility 
(below). 

We choose standard TDGL dynamics driven by gradi- 
ents in the chemical potential, so that the particle current 



W = Min [1, exp (-0AH)] , 



(1) 



-M(<t>)V 



5F 



(3) 



2 



+ 

upward 
triangle 



downward 
triangle 



o 

QO 


o 

QO 


up 


up 


DO 

down 


DO 
Q.„. 

down 




n.n. hops 



n.n.n. hops 




-l 



FIG. 2. 
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normalized asymmetry measure is 
asym = (n up — n down )/(n up + n down ), where the configu- 
rations on top contribute to n up and those below contribute 
to rid 



where M(cj>) is an order parameter dependent mobility. 
A continuity equation is then used to determine the evo- 
lution of the order parameter: 



dt 



-V • J 



E 



dM(4>) 
dz 



(4) 



The choice of a constant mobility M(<fi) — M leads 
to the field dependence dropping out of the dynamics. 
The next simplest choice, M(<fi) — Mn(l — <fi 2 ), the ex- 
act mobility for non-interacting lattice-gases, leads to a 
non-trivial field-dependent DDS coarsening ||JlJ,|l^]. In- 
deed, because the dynamics are deterministic, a semi- 
quantitative understanding can then be reached for the 
linear stability of interfaces and other interfacial proper- 
ties |9|. More generally, we want the mobility to reflect 
the effective coarse-grained mobility. Starting from a 
stochastic model with an applied field, the coarse-grained 
mobility will generally be anisotropic, as can be seen ex- 
plicitly near the critical point S. As a minimal step, we 
allow for different mobilities for currents in the x and z 
directions with 



M x (4>) = (l + m)M (l-tf), 
M z (cf>) = M (l - </> 2 ), 



(5) 



where m describes the mobility enhancement transverse 
to the field direction |Q. In the simulations presented 
here we take a w 2.75, so that the bulk phases are at 
4> = ±0.8. We always set Mq = 1, which fixes the overall 
timescale. The initial conditions <j)(r, 0) follows a Gaus- 
sian distribution around ((f>). 



Asymmetry Measure 

In order to quantitatively compare the models, we need 
a measure of triangular anisotropy. We use the micro- 
scopic measure shown in Fig. |[ That is, we examine all 
squares of four nearest-neighbor sites on our lattice and 
define n up as the number of squares in which the bottom 
two sites are positive but the top two sites have opposite 
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FIG. 3. Configuration snapshots of the Ising DDS with nn 
(nearest neighbor) and nnn (next nearest neighbor) hops. The 
triangles point in opposite directions with respect to the field. 
In both cases, (3E = 0.5 and T = 0.5T C . 



signs. These configurations point "upward" . A similar 
definition is used for ndown- The normalized asymmetry 
measure is then 



asym 



Udo 



< j up 



+ TLdo 



(6) 



so that asym — 1 if all triangles are upward pointing 
and asym = — 1 if all triangles are downward pointing. 
Squares with more or less than three filled sites are not 
counted. The same measure is used for the continuum 
model except we look at four neighboring mesh points 
and count 'full' and 'empty' as cf> > and cf> < 0, respec- 
tively. [In practice an equivalent asymmetry measure for 
the 'empty' phase can be constructed. Similar results are 
obtained.] Our measure is quantitatively different from 
that of Alexander et al Js) where normalization drives 
their asymmetry towards zero as domains grow larger, 
however we obtain the same qualitative sign of the tri- 
angular asymmetry. We prefer our measure since it only 
depends on the shapes of triangles and not their size, at 
least in the coarse-grained formulation. It qualitatively 
agrees with anisotropies seen "by eye" . 



Lattice Gas Results 

By including next nearest neighbor jumps we can flip 
the direction of the triangular domains with respect to 
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FIG. 4. The asymmetry measure as a function of time for 
the Ising model with nn and nnn hops, where (5E — 0.5. The 
measure seems to approach an asymptotic value. 

the field. Fig. | illustrates our results for the Ising model 
with nn and nnn hops. This indicates that the sign of 
the triangular anisotropy is non-universal. We also ob- 
serve that the evolution is more rapid when nnn hops are 
allowed, as probed by domain size. 

The quantitative evolution of the asymmetry as a func- 
tion of time for both nn and nnn hops is shown in Fig. 
The data is averaged over 5 to 10 configurations to reduce 
the noise. The asymmetry starts small and eventually 
saturates to an asymptotic value. We cannot rule out fur- 
ther change [indeed slight decay is evident for T = 0.75TJ 
since there is no known dynamical scaling in the correla- 
tions, i.e. no time-independent scaling function. 

The asymptotic asymmetry vs. T/T c , where T c = T c (0) 
is the critical temperature for zero field, is shown in Fig. 
||. With only nn hops the asymmetry is always positive, 
indicating upward pointing triangles, and increases with 
decreasing temperature. We note that the anisotropy 
measure is continuous across T C (E) (which ranges from 
T c (0) to approximately 1.4T C (0) at E = oo @). At T C (E) 
we would expect surface tensions and particle mobilities 
to be isotropic in a small field, and so we attribute the 
residual anisotropy to the field. [We note that long-range 
correlations can persist into the high-temperature dis- 
ordered phase due to violations in detailed balance in 
driven systems Q .] 

With nnn hops, the triangular asymmetry is small and 
positive near the critical point. This is consistent with 
the previous discussion of nn hops. At lower temper- 
atures the asymmetry turns negative. Since the nnn 
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FIG. 5. The asymptotic asymmetry vs. temperature for 
the Ising model. The nn model asymmetry is always positive, 
while the nnn model asymmetry is positive near the critical 
temperature and at the lowest temperatures. Some field de- 
pendence is also shown. 



hops should lead to a more isotropic particle mobility, 
we might infer that the increasing anisotropy of the sur- 
face tension with decreasing temperature feeds a nega- 
tive triangular asymmetry when not sufficiently 'compen- 
sated' by an anisotropic mobility. Clearly the situation 
is complicated since to fully characterize the 'anisotropy' 
requires the entire function of, e.g., surface tension vs 
interfacial orientation. Without actually having a quan- 
titative measure of the coarse-grained properties of the 
system (surface tension, particle mobility) it is difficult to 
discuss the exact origins of the triangular asymmetry. In- 
deed, this difficulty is a primary motivation to explore the 
coarse-grained picture. Regardless, the particle mobility 
will depend on the microscopic structure, which in turn 
depends on the applied field, even at T = 0. In the limit 
of small applied field E additional induced anisotropics 
should become negligible. We see some indications of 
this through the reduction in the low-temperature posi- 
tive asymmetry regime with nnn hops as the field is re- 
duced, indicating that the regime is induced by the finite 
field. This suggests one possible simplifying tactic: to 
look at the small field limit. Unfortunately this makes the 
timescales for numerical investigation of the driven sys- 
tem inaccessibly large. Analytically, this limit has been 
profitably used by one of us in a coarse-grained analysis 
of surface instabilities in nearly isotropic systems M . 



Coarse-grained Results 

The results for the kinetic Ising model indicates 
that the positive asymmetry measure may be due to 
anisotropy in the mobility. Therefore to flip the trian- 
gles positive in the TDGL model, we allow for differ- 
ent mobilities in the the x and z (field) direction with 
M x /M z = 1 + m. We also varied the bulk coexistence 
value, the field strength E and the initial filling fraction. 
We found that the primary effect came from m and from 
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FIG. 6. Configuration snapshots at t — 1600 and 
t — 12800 of an evolving coarse-grained system with m — 1 
and E = 0.1. The reversal of the asymmetry, while not dra- 
matic, can be seen. The applied field is downward in both 
snapshots. 



E. 

The top snapshot in Fig. || shows that the triangles 
are in the same direction as the n.n. kinetic Ising model 
at early times. This is confirmed by a positive asymme- 
try measurement at these times. However, the bottom 
snapshot in Fig. ^| shows that the asymmetry measure 
becomes negative at late times when the domains are 
very elongated in the field direction. This transient be- 
haviour was robustly present at all nonzero values of m 
and E that we tried. However, the early time "transient" 
regime with positive asymmetry is quite large and can 
be extended indefinitely in the limit of E/m — > 0. This 
is shown in Fig. ^ which shows the asymmetry vs. time 
for critical quenches at fixed m = 1 and varying E. This 
raises the intruiging question of whether the asymme- 
tries seen in the Ising DDS models, where the underlying 
dynamics are slower, might switch at later times. 



Conclusion 
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FIG. 7. The asymmetry as a function of time for fixed 
mobility anisotropy, m = 1 and different fields. The asymme- 
try starts positive but then becomes negative at late times. 
The time for which the asymmetry is positive is extended for 
smaller fields. 



in both stochastic Ising and deterministic coarse-grained 
DDS models. Hence, we see no qualitative differences be- 
tween these approaches with finite fields, and rather see 
great promise in using the strengths of each approach to 
explore DDS phenomenology. 

Questions remain concerning the origins of the trian- 
gular anisotropy within a full anisotropic coarse-grained 
model. This must be understood to intelligently explore 
the parameter space of coarse-grained models. With this 
understanding, we might even profitably turn the tables 
and use the triangular anisotropy as a probe of interfacial 
properties. 

We are not overly concerned with the apparent tran- 
sient nature of the "flipped" anisotropy in the coarse- 
grained model, since even the stochastic Ising model has 
not yet been extensively enough studied to tell if the 
anisotropies hold asymptotically late. However, we will 
now focus our efforts on prolonging the reversal in coarse- 
grained models. This provides a motivation to more fully 
understand the origins of the triangular anisotropy. In 
the process we would like to develop a more intrinsically 
coarse-grained measure of anisotropy that can be used 
equally well in Ising and coarse-grained approaches. We 
expect that anisotropies in the Porod tail of the struc- 
ture factor will be the most robust measure, since they 
directly probe the distribution of interfacial orientations. 
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We have shown that the sign and magnitude of trian- 
gular anisotropies of growing domains are non-universal 
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